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ABSTRACT 

A brief summary is presented of the basic qualitative and quantitative 
aspects of reconnection in its magnetopause setting. First, the basic 
morphological and dynamic features of asymmetric reconnection are examined 
with emphasis on the important role played by the rotational discontinuity 
in these geometries. Second, the structure and other properties of rota- 
tional discontinuities are discussed. Third, the manner in which individual 
particles are energized or de-energized during their interaction with current 
layers in general, and rotational discontinuities in particular, is examined. 
Finally, the question of nor.steady, localized reconnection and its relation 
to flux transfer events is discussed and a qualitative model is proposed to 
describe these phenomena. 


1. INTRODUCTION 


Magnetic field reconnection in a plasma may in principle occur wherever 
the magnetic field exhibits strong shear. In planetary magnetospheres, the 
two principal active sites that have been considered are the magnetopause 
and the magnetotail. In this paper, attention Is focussed on reconnection 
In its magnetopause setting. An overview is presented of the basic features 
and local signatures of the process predicted by existing theory. Important 
new Information concerning magnetopause reconnection has been obtained during 
the last few years, principally as a result of the ISEE mission. These 
results, along with the magnetospheric consequences of magnetopause recon- 
nection will not be dealt with in detail, since they form the topic of 
several papers to follow in this volume. 

2. BASIC MORPHOLOGY AND DYNAMICS 

In the early closed model of the magnetosphere (Johnson, 1960), the 
magnetopause was that surface, usually marked by an electric current sheet, 
which separated the earth's magnetic field from the solar-wind plasma and 
the interplanetary magnetic field (IMF) embedded in it. No interconnection 
between the two fields was included so that the magnetopause was a tangential 
discontinuity. When the concept of magnetic field reconnection was intro- 
duced into magnetospheric physics by Dungey (1961), the meaning of the term 
magnetopause became blurred. Indeed in Dungey 's (1961) drawing of the open 
magnetosphere for purely southward IMF (Fig. 1), the magnetopause seemed 
to be entirely absent, except perhaps near the subsolar point. A geometry 
of this type would arise only if essentially all of the interplanetary magnetic 
flux impinging upon the magnetosphere could be reconnected. 
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The first quantitative analysis of the magnetopause reconnection process 
was carried out by Levy et al. (1964) who concluded that only seme 10-20% 
of the Incident magnetic flux would reconnect while the remainder would be 
carried past the magnetosphere without becoming interconnected with the 
earth's field. In such a situation, the effect of reconnection on the day- 
side magnetopause can be thought of as a small perturbation. In other words, 
the magnetopause remains a well defined current sheet as shown in Fig. 2. 
However, the physical character of this sheet changes drastically: instead 

of being a tangential discontinuity (TO), i.e., a layer with a vanishing 
normal magnetic field component, B n , it is now a rotational discontinuity 
(RD) and has a small but significant B n , amounting to 10-20% of the total field. 

A rotational discontinuity is a large-amplitude Alfvfin wave and the 
interplanetary plasma flows across it with speed equal to the Alfvfin speed 
based on B n . During its passage through the dayside portion of this current 
layer, this plasma is accelerated away from the subsolar point by the 
force, where I^ is The magnetopause current, to form a boundary layer of 
jetting plasma immediately inside the magnetopause. Tailward of the two 
cusp regions, the direction of ^ is reversed so that the plasma is decelerated 
instead (a feature not included in the Levy et al . model). It is in these 
latter regions that mechanical energy is extracted from the solar wind and 
stored as magnetic energy in the geomagnetic tail (e.g.. Swift, 19PJ). How- 
ever, the presence of such deceleration regions is merely a consequence of 
ongoing reconnection somewhere on the subsolar magnetopause. For this reason, 
we shall focus attention on this latter region. 

In the Levy et al. model (Fig. 2), It was assumed that no plasma is 
present in the magnetosphere. As a result, the inner edge of the plasma 
boundary layer inside the magnetopause consists of a narrow slow-mode expansion 
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fan in which the plasma expands to zero pressure and Is accelerated 
further (see Yang and Sonnerup, 1977). In reality the dilute but hot 
magnetospherlc plasma has a substantial pressure so that the expansion 
fan may be absent or even replaced by a slow shock. It Is not clear how 
collision-free slow shocks and expansion fans manifest themselves In the 
narrow magnetopause- boundary layer region. 

The rotational discontinuity Is expected to be present at the magneto- 
pause In all but the most unusual circumstances: strictly antiparallel 
fields and identical or nearly identical plasma and field states on the 
two sides of the magnetopause In which case the usual synmetrlc Petschek 
model containing pairs of slow shocks applies. Even a small Increase In 
density outside the magnetopause over that Inside is sufficient to bring 
the outer slow shocks In that model to their maximum strength (switch-off 
of the tangential field). After that, an RD will appear. If the reconnect- 
ing magnetic fields are not antiparallel, an RD Is always needed. For these 
reasons, the principal magnetic and plasma signatures of the magnetopause 
region away from the separator (reconnection line; X line) should be those 
associated with an RD. 

The separator Is usually assumed to pass through the subsolar point on 
the magnetopause and to be oriented along the net magnetopause current, I, 
as illustrated in Fig. 3a. Figure 3b shows that this orientation does not 
permit reconnection between fields Bj and B 2 , where for an angle as 

such that cosa9>B 1 /B 2 , because the magnetic field component perpendicular to 
the separator must reverse sign for reconnection to be possible (see e.g., 
Sonnerup, 1974). However, Cowley (1976) has argued that the separator orienta- 
tion along j[, while consistent with theory, may not be required by it. A 
hypothesis has also been advanced (Crooker, 1979) to the effect that reconnec- 
tion occurs only at those locations on the magnetopause where B x and B 2 are 


anti parallel or nearly antlparallel (see J. Luhmann, this Volume). The 
merits of this suggestion have yet to be evaluated. 

Direct observational proof of the occurrence of reconnection would 
consist of a measured electric field along a separator since that Is a 
standard definition of reconnection. However, this Is difficult to accom- 
plish, not only because electric-field measurements In the magnetopause 
plasma environment are difficult, but also because the diffusion region, 
i.e., the narrow channel around the separator in which the frozen magnetic 
field condition is violated, has small physical dimensions and therefore 
is difficult to Identify. The difficulty Is compounded by the fact that 
our theoretical understanding of the Important plasma processes In the 
diffusion region Is poor (J.F. Drake, this Volume). In effect, we do not 
know what plasma signatures to look for. A suggestion concerning a magnetic 
signature produced by Hall currents may be found in Sonnerup (1979). 

t 

It is possible to obtain persuasive evidence for reconnection even 
from satellite traversals of the magnetopause away from the reconnection 
site itself, I.e., away from the separator. An electric field E t tangential 
to the magnetopause or a magnetic field component B n normal to it. If present 
over a region of linear dimensions much greater than the magnetopause thick- 
ness, indicates that reconnection is occurring or has occurred in the recent 
past. According to simple steady-state 2D reconnection theory E^ is equal 
to the reconnection rate and B n Is proportional to Both E t and B n are 
difficult to determine: over most parts of the magnetopause, the largest 
electric field component is normal to the magnetopause (see T. Aggson, this j 
Volume) and the largest magnetic field component is tangential to it. Thus 
a reliable determination of E t or B n depends critically upon the knowledge 

of a reliable normal vector. Because of wave motion and other irregularities *] 

» • 

*f 

l 
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In the magnetopause surface, a model normal Is nut useful for this purpose. 
The remaining possibility Is to obtain a normal vector (and B n ) from mini- 
mum-variance analysis of the magnetic field data (Sonnerup and Cahill, 1967). 
Occasionally, reliable normal vectors and B R values are obtained from this 
process and there Is little doubt (e.g., Sonnerup and Ledley, 1979) that 
a few magnetopause crossings have been Identified where B n was significantly 
different from zero. But more often than not, the minimum variance method 
falls to give a reliable normal vector (and B n ) and so the opportunity to 
measure the reconnection rate for a given magnetopause crossing Is limited. 

A third quantity which Is also proportional to the reconnection rate is the 
flow velocity component, v n » across the magnetopause. However, It Is dif- 
ficult to measure and suffers from the same drawback as E^ and 6 n , namely, 
an accurate £ vector Is needed. In addition, a normal flow velocity com- 
ponent Is In Itself not convincing evidence of reconnection since it could 
be the result of rapid diffusion of plasma across a TO. 

A simple check on the presence or absence of reconnection can be obtained 
by examination of the behavior of the tangential plasma velocity across the 
magnetopause. According to the reconnection model, the Interplanetary 
plasma experiences a change in tangential momentum caused by the IxB^ force 
as It crosses the current layer. On the dayside, and with the interplanetary 
magnetic field due essentially south, the net effect is plasma acceleration 
to velocities of about twice the Alfven speed, 2v A , leading to the plasma 
jets shown In Fig. 2, just Inside the magnetopause. These jets have indeed 
been observed (Paschmann et al . , 1979, Sonnerup et al., 1981, Gosling et al., 
1982), and the observed detailed agreement with the theoretical tangential 
momentum change makes it unlikely that these jets were the product of pro- 
cesses other than reconnection. 
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In spite of the fact that the force Is proportional to the recon- 
nection rate via B n , the plasma momentum change Is Independent of that rate. 
The reason Is that the mass flow rate. pv p , across the magnetopause Is also 
proportional to the reconnection rate. Thus, observations of this type 
cannot be used to establish the reconnection rate. 

There are a number of other direct consequences of the reconnection 
field topology which are observable in the region between the Inner and 
outer separatrix surfaces, shown In Fig. 2. These surfaces intersect In 
the diffusion region and they might therefore bear the signature, In the 
form of heat flow, of any electron heating in that region (see J.D. Scudder, 
this Volume). In addition, the outer separatrix might be traced by escaping 
magnetospherlc electrons. 

Between tne outer separatrix and the magnetopause one might expect to 
find magneto spheric and ionospheric Ions that have leaked across the magneto- 
pause and have been given a tangential velocity change of up to 2v A by the 
force (Scholer et al . , 1981; Sonnerup et al., 1981), as well as magneto- 
sheath ions that have been reflected In the magnetopause and have been simi- 
larly influenced by the I*B, force (Sonnerup et al., 1981). Between the 

« I 

magnetopause and the inner separatrix, one expects to find, not only the 
jetting plasma boundary layer, but also any magnetospheric ions that have 
not leaked across the magnetopause but have been reflected against it (Scholer 
and Ipavich, 1983). Detailed consideration of these various signatures 
may be found in G. Paschmann's article (this Volume). Finally, it is noted 
that the inner separatrix comprises the last set of closed field lines in 
the magnetosphere. This surface and the region just outside it project into 
the high-latitude the ionosphere so that many of the consequences of magneto- 
pause reconnection should be observable there (see P. Reiff, this Volume). 
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3. ROTATIONAL DISCONTINUITY 

The basic properties of an RD are: (1) It has a nonvanishing magnetic 
field component, 8 n , normal to the layer; (11) the plasma flows towards 
and away from the discontinuity with normal speed, v n , equal to the Alfven 
speed based on B n and corrected for nonlsotroplc pressure; (111) the tan- 
gential magnetic field can change direction by an arbitrary anqle ae across 
the sheet; Indeed, It Is the only discontinuity having this capability. 

The fact that the flow speeds Into and out of the discontinuity are equal 
to the corresponding Alfven-wave speeds guarantees that no wave steepening 
or broadening of the usual type Is present In a uniform medium. 

Jump Conditions 

In a collision-free plasma with nonlsotroplc pressure the jump con- 
ditions across an RD are (Hudson, 1970; 1971; 1973): 


Mp(l-a)} * 0 

(Mass Cons.) 

0) 

MP i +B 2 /2u 0 ) - 0 

(Normal Mom. Cons.) 

(2) 

A{ vV * 0 

(Tang. Mom. Cons.) 

(3) 

A{(a+%)B 2 /u Q p+5Pj/2p} * -q 

(Energy Con».) 

U) 


In these equations, the symbols, p, p, B, and v^ have their usual meaning 
and the delta bracket is defined by a{K) s K 2 -K , where the subscripts 1 
and 2 denote conditions upstream and downstream of the discontinuity. Also, 
q Is the amount of energy per unit mass which leaves the system via heat 
conduction or radiation. The pressure anisotropy, o, is defined by 


* 5 (P| -Pj^Wq/B 2 


(5) 


8 


Finally* y At = B t C ( 1 -a }/y Q p is the tangential Alfvftn speed. In the frame 
of the RD, the normal plasma flow speeds towards and away from the layer 

are v XJ s v Axi » B x ^ 1 “ 0l ^ w o p i^ and V X 2 " v Ax 2 * res P ectivel y* I n the fire- 
hose limit, °i * <*2 * 1» the P lasma ceases to flow across the discontinuity. 

From the above formulas with q»0 and aj'a^O, the well known MHD results 
(e.g., Landau and Lifshitz, 1960) are recovered: p, p, and B remain un- 

changed across the RD while the direction of B_ t can change by an arbitrary 
angle, this change being accompanied by a corresponding change in v t as 
indicated by Eq. (3). If the tangential field rotates by an angle A8, 

Eq. (3) yields |av^| * 2jv At | s i n( A0/2) so that the maximum value of |av^| 
is 2v At , a result already quoted in Section 2. 

If q=0 and a 1 «a 2 ?3/4 t p, p, and B again remain constant but the AlfvGn 
speed now contains the correction factor (l-a)* 5 . This situation arises if 
the double adiabatic relations Pj/pB = const, and p fl B 2 /p 3 = const, hold. 

If only the first of these conditions, representing the conservation of the 
magnetic moment of a particle, is valid, a curious situation arises (Hudson, 
1973). One possible root is a 2 =a 1 in which case p, p, and B all remain 
unchanged. But there also exist parameter ranges in which a second root 
a 2 ^. x occurs in which case p, p, and B all change across the RD. In discon- 
tinuities of this type the entropy usually also changes: for q=0 only solu- 
tions that bring about an entropy increase, A{s}>0, are physically acceptable. 
If neither of the adiabatic relations hold and/or if q^O, even more general 
behavior appears possible. 

In summary, contrary to the case of shocks, the jump conditions across 
a rotational discontinuity do not uniquely specify the change of state of 
the plasma and field. As long as q=0, one possibility is always that the 
plasma state *-.id field magnitude remain unchanged but other possibilities 
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exist as well. These have been explored to some extent by Hudson (1971, 

1973). However, recent experimental evidence (6. Paschmann, this Volume) 
indicates that substantial ion heat flow away from the magnetopause RD 
occurs occasionally so that the assumption q*0 Is not always a good one. 

Structure 

The structure and thickness of the rotational discontinuity must be 
such that the changes implied by the jump conditions are achieved. For 
example, one must presume that the two conditions q*0 and Ms}=0 together 
imply a sufficient thickness of the layer so as to permit laminar, non- 
dissipative behavior of the plasma. It also seems likely that limitations 
derived from the RD structure may eliminate certain downstream states that 
are allowed by the jump conditions. 

In dissipationless MHD, the RD structure is simple: the plasma state 

and the magnetic field magnitude remain constant while the angle, e, of the 
tangential field changes by the desired amount, ae. The function e(x), x 
being the normal coordinate, can be specified arbitrarily. If a small amount 
of dissipation is present, the width of the layer increases gradually. 

Our theoretical knowledge of the structure of rotational discontinuities 
in a collision-free plasma is not extensive and is limited to laminar struc- 
tures, as discussed below. However, observations (e.g., Sonnerup and ledley, 
1979; Berchem and Russell, 1982; Paschmann et al., 1979) suggest that in 
reality the structure may be turbulent. 

In analyzing the structure of current layers, it is convenient to use 
tie so-called de Hoffmann-Teller (dHT) frame in which the external electric 
field vanishes (such a frame can be found for all plane one-dimensional layers 
except perpendicular shocks). In this frame, the Inflow and outflow velocities 

4 

at an RD are field aligned and equal to the AlfvSn speeds on the two sides. 
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The following basic statements can be made: 

(1) The structure of a rotational discontinuity involves an electric 
field E n (x) normal to the layer and an associated potential barrier $(x). 

The first-order orbit theory model examined by Su and Sonnerup (1968) as 
well as the double-adiabatic model by Lee and Kan (1982) contains such an 
electric field. But it can be shown in a general way (Sonnerup and Wang, 

1984) that, regardless of model, an RD must always contain a potential 
barrier. Su and Sonnerup recognized the possibility of trapping particles 
electrostatically in the layer but they did not satisfy charge neutrality 
(or Poisson's equation). Lee and Kan were the first to understand the 
crucial role played by trapped electrons in this regard and theirs is the 
only complete model in print to date. It predicts E n values of the order 
of one millivolt per meter for typical magnetopause conditions. However, 
recently Wang and Sonnerup (1983a,b) have developed a model in which large 
deviations from charge neutrality occur in narrow regions (electrostatic 
"shocks") of width equal to a few Debye lengths. The electric fields in 
these narrow regions are of the order of volts per meter. 

It should be noted that both E n (x) and e(x) change when one transforms 
from the dHT frame to some other frame of reference (the lab frame) moving 
parallel to the discontinuity. This must be kept in mind when comparing 
electric field measurements to theory. 

(2) For an RD in a plasma consisting of electrons with isotropic pressure 
tensor and one specie of ions of charge q, the tangential magnetic field, 
B t =By+1B z , where i 2 =-l, can be shown to obey the equation (Sonnerup and 

Wang, 1984) 


dB f /T-aT 
35T* « t Tq7 



( 6 ) 


11 


where X^ = ( m i/w 0 n ii q2 ^ ls the ^ on ^ nertial length, is the ion number 
density, and is the pressure anisotropy factor, upstream of the RD. Also 
P ixt 2 P 1xy + ^ P ixz the tangential ion stress. Note that because 
^ (1-ai) * n^(l _a 2 ) all upstream conditions (subscript 1) can be replaced 
by downstream ones (subscript 2). It is easy to shew that for the upstream 
or downstream state, P^ xt * B x B t a/y Q so that dB^/dx = 0 there. 

In the absence of viscous stresses, the ion pressure tensor P. is diagonal 
in a coordinate system with one axis along £. In that case P^ xt is propor- 
tional to B fc so that the ratio p -j x t/ B t ls P urel y real. Equation (6) is 
then of the form dB^/dx = ig(x)B^., where g(x) is a real function of x, indi- 
cating that the field magnitude |B t | is constant. 

In magnetopause rotational discontinuities the magnetic field magnitude 
often has a minimum in the center of the layer. The above results indicate 
that such an effect is likely to be caused by viscous stresses in the layer 
(although nonisotropy of the electron stress tensor and/or the presence of 
more than one ion specie may be contributing factors as well). Viscous 
stresses are expected to be important only in thin layers so that strong 
deviations from |B t |=const. should be an indication that the magnetopause 
width is small. 

(3) Equation (6) also forms a suitable basis for discussing the sense 
of polarization of the RD. The right-hand (electron) and left-hand (ion) 
polarizations are obtained when the rectangular bracket is positive and 
negative, respectively. For example, if the ions are cold so that P 1xt =0, 
<^“0, then the ion polarization is obtained for a positive potential barrier, 
<f>(x)>0, since such a barrier will lead to n^>nj . Similarly, $(x)<0 yields 
the electron polarization. For hot ions, the situation is less well under- 
stood. For example, in the double-adiabatic description of the ions employed 


12 


by Lee and Kan (1982) the square bracket reduces to [l-n^/n^] but for <t>(x)>0 
(the only case dealt with by Lee and Kan) the density ratio n^/n^ is found 
to be less than unity except for cold or almost cold ions, the result being 
the electron polarization. 

Sonnerup and Ledley (1979) have argued that only the electron polariza- 
tion should occur for layers that are sufficiently thin so that only elec- 
trons but not ions are capable of moving across the layer by sliding along 
the magnetic field lines so as to provide the required field-aligned current 
distribution in the layer. This situation should arise when the layer 
thickness is comparable to, or less than, the ion gyroradius. On the 
basis of Explorer 12 and 0G05 data, these authors also state that both 
polarizations have been observed but with a preference for the electron 
sense. On the other hand, Berchem and Russell (1982) (see also Russell, 
this Volume) have found from ISEE data that the basic polarization rule 
is for the field rotation angle A6 in the layer to obey the inequality a9<tt. 
However, these authors did not distinguish between rotational and tangential 
discontinuities so that the distribution among the former between the ion 
and electron polarization in their data set remains unknown. However, 
examination of published Explorer 12 and 0G05 RD's basically supports the 
Berchem and Russell rule. Only one crossing has been published in which 
A6>tt (Sonnerup and Ledley, 1974). This crossing had the electron polariza- 
tion. 

Further support for the rule A6<tt comes from a recent computer simula- 
tion of RD's by Swift and Lee (1983) in which cases with Ae>Tr were observed 
to be unstable. The reason for this effect is not clear. 
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4. PARTICLE ORBITS 

Motion In de Hoffmann-Tel ler Frame 

In this section we review certain general results concerning particle 
orbits in one-dimensional current sheets with 8 x *B n fO. It is advantageous 
to study the orbits in the dHT frame because in that frame particles have 
simple helical orbits as they travel towards and away from the current 
sheet. In the sheet itself their motion may be extremely complicated but 
it renal ns constrained by the conservation of total energy and the two 
generalized tangential momenta: 


mv 2 /2+q$(x) « const. 

(7) 

mv -qB z+qA (x) = const, 
y * y 

(8) 

mv z +qB x y+qA 2 (x) * const. 

(9) 


where = -dA z /dx and B z 5 dA^/dx are the field components tangential to 
the layer. 

Cowley (1978) has shown that Eqs. (7), (8), and (9) constrain the 
particle orbit to lie inside a surface in space, the intersection of which 
with any plane parallel to the yz plane is circular with radius R x * mv/qB x - 
As one moves from one plane parallel to the yz plane to another, R x changes 
as the particle speed v changes in response to <j>(x). The center of the 
circle moves in such a way that it remains located on one and the same 
field line. 

It is instructive to arrive at these results by use of the diagrams in 
Fig. 4. As shown in part (a) of that figure, at a given fixed x value, the 
position of a particle in velocity space is constrained to lie on a sphere 
of radius v where v = v /l~2q$(xj/m, being the particle speed upstream of 


the current layer where $=0. On this sphere the circles which represent 
constant pitch angles, a, with respect to the local magnetic field are shown 
for convenience. If the geometry shown in Fig. 4a is projected onto the v 
v 2 plane, the result is the circular disk v 2 + v 2 <v 2 shown in Fig. 4b. The 
circles of constant pitch angle now appear as ellipses inside the disk. The 
cases a*0, corresponding to particle motion along B towards increasing x, 
and a=ir, corresponding to motion along B towards decreasing x, appear as two 
points labeled, 0 and n, respectively. For a fixed x value, Eqs. (8) and 
{9) provide a linear transformation from the tangential velocity map in 
Fig. 4b to a plane in configuration space parallel to the yz plane and 
located at the assumed x value. This transformation consists of a clockwise 
rotation of the disk and a scaling by the factor m/qB x , as shown in Fig. 4c. 

The terms A^(x) and A z (x) in the transformation lead to a shift of the 
disk location with changing x value such that the center remains on one 
and the same field line. To see this, we differentiate Eqs. (8) and (9) 
to obtain 


B x dz = (dAy/dx)dx = B z dx 

(10) 

B x dy * (-dA z /dx)dx * B^dx 

(11) 


which are the differential equations for a field line. 

In the uniform field region on either side of the current layer, the 
surface traced by the disk (Fig. 4d) is a cylinder of elliptical cross 
section as shown in Fig. 4e for the upstream side, say, of the layer. In 
the figure, the symbol f denotes the angle between the incident magnetic 
field, assumed to lie in the xy plane, and the surface of the discontinuity. 
Inside the ellipse are circles with centers at z * (mv/qB x )cosJ^ cosa, each 
labeled by a corresponding pitch angle, a. For 0<a<ir/2 the particle moves 
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In a helical orbit toward the layer while for ir/2<a<ir it moves away from it. 
The two foci of the ellipse correspond to a*0 and a*ir. For pitch angles 
near 0 and *, the circles do not touch the ellipse: they correspond to 

helices in which v x does not change sign. For circles that touch the 
ellipse. v x changes sign at the points of osculation. 

As a simple illustration of Fig. 4e, consider a particle incident upon 
the current sheet with pitch angle a x (<tt/2) which is reflected and leaves 
with pitch angle a 2 (>ir/2). Figure 4e shows that during the interaction 
with the current sheet, the guiding center of the particle is displaced 
by the distance 

az * jjjjp cosf^ (cosaj-cosa 2 ) (12) 

i.e., purely in the z direction. The maximum displacement occurs when 
a x a 0, a 2 B Tr, in which case the orbit moves from one focus of the ellipse to 
the other (a result that may be compared to the approximation used by T.W. 
Speiser, this Volume). 

A second illustration is provided in Fig. 5 which shows a view along 
the normal, i.e., along the positive x direction of a rotational discontinuity 
with A0*135°. A field line bends, as shown in the figure, as it passes 
through the current layer. Disks are shown at the entrance to (#1), in 
the middle of (#l*s), and at the exit from (#2) the layer. Note that the 
constant pitch angle pattern in the disk rotates as one progresses through 
the current layer. Disks #1 and #2 are shown with the same diameter, corres- 
ponding to an assumed value <f> a 0 of the electric potential on both sides of 
the layer. Disk #1^ has a smaller diameter corresponding to the assumed 
presence of a potential barrier #(x). As an example, a particle with zero 
pitch angle enters the layer at the point a=0 in disk #1. If the magnetic 
moment u is preserved, it leaves at the point a a 0 in disk #2. Note that 
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because of Inertia drift it does not follow the field line through the 
entry point which leaves the layer at the point labeled P. If the particle 
were reflected rather tha.i transmitted it would leave at the point a*ir in 
disk #1 instead (assuming constant y). In reality, the magnetic moment is 
not necessarily preserved and numerical orbit calculations are needed to 
obtain the exit pitch angle. 

Particle Energization 

Equation (3) describes the tangential momentum change of the plasma 
as it crosses an RD. Depending on circumstances, this change will lead 
either to energization or to deenergization of the plasma. These effects 
can be examined by transforming velocities from the dHT frame to the 
"laboratory" frame. As illustrated in Fig. 6, this transformation con- 
sists of adding a constant velocity ^ parallel to the sheet. Energization 
of transmitted and reflected particles is illustrated in Figs. 6b and 6c 
while an example of deenergization is shown in Fig. 6d. For simplicity, 
only the case where Bj, B^, and v^ all lie in the xy plane is shown. 

In the laboratory frame a constant tangential electric field, 

E+ = -yxB, will be present, the direction of which is shown in Fig. 6. 

w ■““Tl 

This is the reconnection electric field, and it now becomes clear that 
E*. and B n are indeed proportional. At the magnetopause is a remnant 

of the interplanetary electric field (IEF) and is directed from dawn to 
dusk for southward interplanetary magnetic field (IMF). Note that £ t *I>0 
for energization and £ t *I<0 for deenergization. The former case occurs 
on the frontside magnetopause where the magnetopause current, I_, ts directed 
from dawn to dusk for southward IMF. The latter case occurs on portions of 
the magnetopause tail ward of the cusp where _I is directed from dusk to dawn. 
The suggestion has been made (Crooker, 1979) that it may also occur on the 
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frontside magnetopause as a consequence of cusp reconnection. But since 
I_ Is directed from dawn to dusk on the frontside, E^ would then have to be 
directed from dusk to dawn, l.e.. It would have to oppose the IEF. This 
appears unlikely. 

The energization or deenerglzation of Individual particles occurs as a 
result of a displacement of the particle guiding center along E t - As an 
example, consider a particle reflected at the current layer. It Is easy 
to snow from the velocity triangles In Fig. 6c that the kinetic energy 
increase is 

A e * mvv o cos^(cosa 1 -cosa 2 ) s qE^Az (13) 

(Note that a 2 >ir/2) where az is the guiding center displacement. Since 
v = E t /B * E t /B we can solve this expression for az, the result being 
exactly the displacement given by Eq. (12). In other words, orbit considera- 
tions based on the conservation of energy and generalized tangential momenta 
lead to exactly the same result as the frame transformation method described 
in Fig. 6. This Is not surprising but It serves to emphasize that the MHO 
reconnection model which predicts plasma acceleration by the IxB^ force is 
internally consistent: the reconnection electric field is precisely the 

field needed to energize the plasma particles by the requisite amount. Thus 
experimental checks of the tangential momentum balance automatically provide 
a check of the electromechanical part of the energy balance. However, the 
remaining, thermal part is in general not negligible (see G. Paschmann, 
this Volume). 

Equation (12) indicates that the particle displacement Az along E t 
becomes very large as B x *B n approaches zero. However, at the same time E^ 
itself becomes small, the result being that the energy gain a e remains the 


same. Therefore, as stated earlier the plasma energization is independent 
of the reconnection rate. It is also evident from Eq. (12) that the ener- 
gization of electrons is negligible compared to that of the ions. 

The general formula for the energization of a particle transmitted 
through an RD, as illustrated in Fig. 5, is not much more difficult to 
work out than the simple case of reflected particles discussed above. 

For example, a particle which crosses the layer while retaining a*0 will 
undergo an effective tangential displacement from point P to point 
a*0 in disk §2 so that its energy increase will be Ae * The 

reason that the distance from the point o*0 in disk #1 to the point P 
is not included in Ar t Is that these two points have the same electric 

potential: the total transformation electric field E„ * -v xB has no 

— o — o — 

component along B and the two points are located on one and the same 
field line. However, as mentioned already in Section 3, E^ does have a 
substantial component along the normal direction. 
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5. NONSTEADY LOCALIZED RECONNECTION 

The preceding sections have dealt with reconnection as a steady-state 
process. However, there Is mounting evidence (Russell and Elphlc, 1978; 
Cowley, 1982) that a nonsteady patchy version of reconnection, referred to 
as flux transfer events (FTE's), may be the dominant reconnection mode at 
the magnetopause. The FTE Is envisaged as a pair of flux tubes, one In 
each hemisphere, each passing through a "hole" In the magnetopause and 
connecting to the earth and each being accelerated along It away from the 
subsolar region by the magnetic tension associated with the sharp kink 
In the tube at the hole (see Fig. 7a). Such tubes seem to have typical 
cross-sectional dimensions of the order of 1 (Implying a hole dimension 
of a few R £ since B n <B) but In reality there may be a continuum of sizes 
ranging from dimensions comparable to the Ion gyroradlus up to the scale 
of quasisteady reconnection. Detailed discussion of the morphology and 
statistics of FTE's may be found elsewhere In this Voljr.it;. Here a few 
general remarks are offered concerning these structures and their relation 
to reconnection. 

(1) As illustrated In Fig. 7b, the pair of holes in the magnetopause 
associated with FTE's must be the result of reconnection, limited to a 
narrow longitude Interval, which started near the equatorial region and 
then ceased almost immediately before a steady state could be established, 
the method of cessation being conversion of the reconnection line Into two 
passive or almost passive X lines and an 0 line. 

(2) If quasi-steady reconnection takes place In a limited longitude 
segment then, for aefTr, there are regions on the magnetopause surface where 
one would observe open field lines but not the plasma jetting discussed 
previously. Signatures of this type should be looked for In the data and 


their relation to FTE's should be examined In detail. 

(3) Whenever the IMF Is due south, the magnetopause situation would 
seem to be an Ideal one for quasi steady reconnection to occur at Its maxi- 
mum permitted rate. Yet, this Is apparently not what happens. Even with 
a southward IMF, quasisteady reconnection Is observed only occasionally. 
Indicating that the process occurs sporadically and/or that It Is limited 
to narrow longitude segments. If the currently popular interpretation of 
FTE signatures Is correct, then what Is observed most of the time is there- 
fore a patchy time-dependent version of reconnection. This set of circum- 
stances Indicates the existence of a threshold (other than a southward IMF) 
for the onset of the process. Either by design or by accident the magneto- 
pause plasma state seems to hover around this threshold. 

(4) The detailed nature of the threshold Is not understood but the 
following scenario Illustrates how FTE's might be generated. Assume that 
a flux tube containing interplanetary magnetic field lines drapes over the 
subsolar magnetopause and gets hung up there, perhaps In a preexisting 
Indentation. The plasma in this tube will then escape by flowing tallward 
along the lines of force, the result being a lowering of plasma density, n, 

6 value (In particular ^ )» and Alfven Mach number, M^, of the flow. It 
may be argued that each of these factors Is conducive to the onset of re- 
connection between the field lines in the flux tube and the geomagnetic field. 
As soon as reconnection has been initiated, two developments occur: (1) a 

deepening of the Indentation in the magnetopause associated with the flux 
tube must take place as reconnection erodes flux from the magnetosphere; 

(11) the region occupied by the original flux tubs gets replenished with 
fresh solar-wind plasma in which n, b, and M^ return to their original values. 
Via the threshold, the latter effect may lead to the cessation of reconnection. 
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The former effect creates a suitable trap for new Interplanetary magnetic 
flux tubes. These Indentations and reconnection regions may also get 
swept away over the flanks of the magnetopause so tnat FTE’s should be 
observable over the entire frontside magnetopause and not just near local 
noon. And on occasion, the 0-type null line created when reconnection 
ceases (see Fig. 7b) may be swept over one or the other of the poles so 
that FTE signatures characteristic of the southern hemisphere may occasion- 
ally oe observed In the northern hemisphere and vice versa. 

It also seems reasonable to assume that for suitable plasma conditions 
a quasisteady reconnection configuration may be established which either 
remains limited to a narrow longitude segment or spreads over a substantial 
part of the frontside magnetopause. 

The above scenario Is by no means unique, but It may serve as a 
useful guide for future theoretical and observational studies of FTE's. 

In particular. It Illustrates the Importance of developing a better under- 
standing of the threshold conditions for onset of magnetopause reconnect lent. 
It seems clear that these thresholds must be associated with local conditions 
at the reconnection site rather than with global boundary conditions. A 
complete theory is needed for the onset and evolution of the tearing mode 
in a collision-free plasma when the magnetic field exhibits strong shear 
as it usually does in the magnetopause current layer. 
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FIGURE CAPTIONS 

Fig. 1 The first reconnection model of the magnetosphere (Dungey, 1S61). 

The reconnection electric field E t , which leads to the flow pattern 
shown, is comparable in size to the interplanetary electric field 
so that most of the Incident magnetic flux becomes Interconnected 
with the earth's field. 

Fig. 2 Magnetopause reconnection model due to Levy et al . (1964) in which 
only a small portion of the incident magnetic flux is reconnected. 
The principal geometrical features of this model are: the separator 

(X); the magnetopause current layer in the form of a rotational dis- 
continuity (RD); a thin high-velocity plasma boundary layer (BL) 
immediately inside the magnetopause and, at its innermost edge, a 
slow-mode expansion fan (SEF; shown shaded). The inner and outer 
separatrix surfaces are marked by (IS) and (OS), the reconnection 
electric field by E^, and the magnetopause current by I. 

Fig. 3 (a) The separator (X line) is aligned with the net magnetopause 

current I. (b) For cosa 6>B 1 /B 2 no reconnection is possible. 

Fig. 4 (a) In the de Hoffmann-Teller frame, and at a fixed x value, a 

particle is located on a sphere of radius v in velocity space. 
Circles of constant pitch angle on this sphere are shown. 

(b) Projection of the sphere and the constant-pitch-angle 
circles on the v^v 2 plane, (c) The generalized tangential 
momenta provide a linear transformation from the v^v z plane to 
the yz plane, (d) Surface in space traced by the disk (c) as 
x is changed, (e) Elliptic cross section of the surface (d) on 
left-hand side of current layer. 
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Fig. 5 View of the magnetopause from the sun. The interplanetary field 

connects with the magnetospheric field B 2 across the magnetopause. 
Cowley's (1978) disks are shown outside (#1), in the middle of (#l*s), 
and inside (#2) the magnetopause. 

Fig. 6 Velocity change of a particle as it interacts with the magnetopause 
current layer (shaded), (a) In the de Hoffmann-Teller frame, a 
transmitted particle has velocities v : * vcosa 1 and v 2 * vcosa 2 
before and after the interaction, (b) In the laboratory frame 
a transformation velocity ^ has been added and the particle now 
has velocities vj and v 2 instead where v 2 »v]. (c) A reflected 

particle is similarly energized, (d) De-energization is achieved 
by reversing the direction of y 0 . The magnetopause current 2 
and the electric field E. t * are P aral ^ el for energization, 

antiparallel for de-energization. 

Fig. 7 (a) View from the sun of a flux transfer event (FTE) on the mag- 

netopause. (b) Side view of the magnetopause (shaded) during the 
development of an FTE. At t=0 the frontside magnetopause is closed, 
with an 0-type magnetic null line at the subsolar point. Reconnection 
starts, splitting the 0-line into two 0-lines and an X line, the 
latter remaining at the subsolar point. At t*l all the magnetic 
flux comprising the original magnetopause has been reconnected. 

At t=2 some interplanetary magnetic field, B 1# has become connected 
to the earth's field, B 2 . At this time, the X line is converted to 
an 0-line and two X lines, the former remaining at the subsolar point. 
A small amount of reconnection occurs at the two X lines, the result 
being the formation of a closed magnetopause near the subsolar point 
at t=3, along with two separate holes in the magnetopause. The situ- 
ation at t=3 corresponds to the section A-A in part (a) of the figure. 
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ABSTRACT 


According to the usual magnetopause reconnection model, plasma flows 
across the magnetopause current sheet, which is a rotational discontinuity 
with a nonzero normal magnetic field component B p , from the magnetosheath 
into the magnetcspheric boundary layer. As the plasma crosses the sheet, 
which has net current I, it Is accelerated by the IxB„ force and flows toward 
the poles with speeds up to twice the AlfvSn speed. In the past, a number * 

of ISEE magnetopause crossings in which such plasma jetting was observed 
have been Interpreted In terms of reconnection, on the basis of a quantita- j 

tive check of the tangential stress balance across the magnetopause. In this | 

paper, we examine the total energy balance instead for two events {August 
9 and September 8, 1978), with the objective of obtaining an additional check 
on the interpretation In terms of reconnection. The results are: (1) To I 

within experimental uncertainties, the plasma and magnetic field data are * 

consistent with reconnects.:; (11) An enthalpy increase comparable to the 
kinetic energy increase occurs In the magnetopause. Thus substantial dissi- 
pation Is present In the rotational discontinuity; ( i i 1 ) During the September 
8 event, an Ion heat flow associated with a beam of reflected magnetosheath 
particles carried away some 20% of the total converted electromagnetic energy; 

(iv) The energy balance provides a persuasive check of the inferred location 1 

of the ISEE spacecraft as being south of the reconnection line on August 9, 
and north of It on September 8. 



1 . INTRODUCTION 


Three recent articles (Paper I: Paschmann, et al. t 1979; Paper II: 
Sonn^rup, et al., 1981; Paper III: Gosling, et a«., 1982) have contained 
Information on a total of 12 Instances of large plasma flow velocities 
near the magnetopause, observed by the spacecraft ISEE 1 and 2. These 
observations were Interpreted In terms of the magnetic field reconnection 
process, one of the most striking directly observable effects of which Is 
the generation of high-speed plasma jets directed away from the reconnection 
site, l.e., away from the so-called X line. It was shown In Papers I-III 
that, within experimental uncertainties, the observed plasma jet velocities 
In several of the events was consistent with the reconnection model . In 
this model, described In detail In Paper II, the jets consist of magneto- 
sheath plasma that has crossed the magnetopause and has been accelerated 
by the I x ^ force there. Here ^ and are the current 1r, and the mag- 
netic field normal to, the magnetopause, respectively, and the magnetopause 
away from the reconnection line Is a rotational discontinuity. Additional 
support for the Interpretation of these events In terms of reconnection has 
been obtained by examination of energetic particle anisotropies (Scholer, et 
al., 1981; Paper II) and particle reflection at the magnetopause (Paper II; 
Scholer and Ipavich, 1983.) 

On the other hand, the Interpretation cf one event (September 8, 1978) 

In terms of reconnection has been questioned by Eastman and Frank (1982). 
Further detailed observational Information concerning this event may be found 
In Scholer, et al . (1982) and a possible resolution of the observational 
dilemma discussed by Eastman and Frank has been proposed by Daly and Fritz 
(1982). 
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In addition to the tangential momentum balance, which was studied In 
Papers I-III, It Is also desirable to check the energy budget during these 
plasma acceleration events. Indeed, Helkklla's (1975) challenge of the 
magnetopause reconnection model was based on an energy, not a momentum 
argument. The question Is whether the electromagnetic energy, £•£, con- 
verted per unit area of the magnetopause during reconnection can be accounted 
for In terms of the increase In kinetic energy of the magnetosheath plasma 
as It flows across the magnetopause Into the jets. In terms of an Increase 
in the enthalpy of this plasma, and in terms of heat flow away from the mag- 
netopause. It Is worth noting that except for heat flow, ohmic and viscous 
dissipation, the energy conservation law describes reversible thermo- and 
electromechanical conversion and can be derived directly from the momentum 
conservation law and the double adiabatic relations. Thus the momentum 
balance check In Papers I-III automatically provides a partial check of 
the energy budget. 

It Is the purpose of this paper to examine the full energy balance 
at the magnetopause In order to see if a further check can be obtained of 
the Interpretation of the plasma acceleration events In terms of the recon- 
nection model. In Section 2, we review the energy balance equation and 
cast It In a form suitable for application to the ISEE 1 and 2 plasma and 
magnetic field observations. The former are obtained from the MPE/LANL 
fast plasma analyzer, the latter from the UCLA magnetometer onboard the 
spacecraft. In Section 3 we examine two acceleration events from Paper II 
(September 8, 1978; August 9, 1 978) in detail, and In Section 4 we discuss 
the results and present our conclusions. 
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2. ENERGY EQUATION 

The steady-state energy balance In a one-dimensional plasma discontinuity 
may be written In the form 

o (y_*n_)A( I 5V 2 +U} + A{y_*P*£} + A{£»£) + a{ (£*B/u 0 )*£} * 0 (1) 

Here p and v_ are the plasma density and velocity, respectively, and £ is the 
heat flow vector. Also, P Is the full electron and Ion pressure tensor, U Is 
the Internal energy per unit mass, and £xB/u 0 Is the Poyntlng vector. The 
symbol a{ } Indicates the change In the quantity In braces from one position 
of observation to another along the vector £ normal to the layer. In the 
present application, £ is the magnetopause normal and Is taken to point 
outward from the earth. Thus the (constant) mass flux across the layer, 
p(v;£). Is negative. In the paper we shall employ a right-handed cartesian 
coordinate system. (x,y,z) with the x axis along £. The two positions of 
observation will be a fixed reference point In the magnetosheath near the 
magnetopause, denoted by the subscript 1, and a movable point In the magneto- 
pause current layer or in the boundary layer Inside It, denoted by the sub- 
script 2. With this notation, we have, for example, A{(J> » *J -U . 

The first term In Eq. (1) represents the change in kinetic and internal 
energy as the plasma moves inward through the magnetopause current layer. 

The Internal energy may be written 

U » (TrP)/2p • /2p+p ± /p (2) 

where p^ and are the pressures parallel and perpendicular to the magnetic 
field, respectively. Thus the incernal energy, l.e., the sum (Pjj/2-iP^), can 
be obtained directly as the trace of a measured pressure tensor but to ‘alculate 
Pjl and p^ Individually, certain assumptions are needed concerning the viscous 
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part of the tensor. 

The second term in Eq. (1 ) Is the net rate of work done on the slab of 
plasma located between points 1 and 2 by the various stress components. It 
Incorporates both the flow work, i.e., the work produced by a diagonal stress 
tensor with components (Pj| ,P^,P^), and the work done by the viscous stresses 
acting on the surfaces of the slab. If we denote the latter by p(v»ji)A{y} 
we may then demonstrate that 

A{y_*P*Jl} * p(v.n)A{p i /p} + (B/ji)A{a(v/Bj/y 0 } + p(y/n)A{J?} (3) 

where a is the anisotropy factor, i.e., a = (p^ -p^)u Q /B 2 . The rate of work 
done by viscous stresses at station 1 Is zero, i.e, * 0, since this point 

is located outside the region of velocity shear. This is also true for 
station 2, w „ = 0, when this point is located in a uniform plasma region 
inside the magnetopause. In that case, no viscous term remains. But if 
station 2 is located in a region of high velocity shear in the magnetopause 
a substantial viscous contribution could be present. It is difficult to 
estimate the magnitude of this contribution but it should be possible to 
determine its sign by use of the regular expression for the viscous stress 
tensor in a newtonian fluid with viscosity y. The result is 

p(v-n)V 2 - -n 2 [^- (v2/2+v|/2+2v|/3)] (4) 

2 

During these events, the tangential velocity components v and v 2 in the 
magnetopause are much larger than v x and they decrease with increasing x. In 
other words, the right hand side of Eq. (4) is positive. Since (v/n)<0, 
we therefore conclude that jy is negative in the magnetopause. 

It is convenient to write the third term in Eq. (1) in the form 


A{£*n} = p(v_*n,)Q 


(5) 
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where Q represents the net Ion and electron heat flow per unit mass away 
from the plasma slab located between points 1 and 2. It Is difficult to 
predict the magnitude of Q theoretically but since the magnetopause and 
boundary layer plasma appears to be heated one would expect a net heat flow 
away from the entire region, l.e., Q>0. Furthermore It can be argued that 
Q should be smaller the closer point 2 Is to point 1. We shall show from 
the observations that this term sometimes makes an important contribution 
to the energy balance. 

Finally, the fourth term in Eq. (1) is the net transport of electro- 
magnetic energy Into the region between stations 1 and 2. Since the tan- 
gential component of the electric field remains constant in the layer and 
since the change In tangential magnetic field Is proportional to the total 
electric current 1^ In the slab between stations 1 and 2, this term may 

also be written in the familiar form E*I . 
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If we write 01m' s law as 


£ + v x b = 


( 6 ) 


where ^ is an effective resistivity tensor and the current density, we may 
transform the fourth term in Eq. (1) to the following form: 


M(ExB)‘n/u 0 } * A{-[(vxB)xB]*n/u Q } + A{[(^*i)xB]*n/u 0 } 


(7) 


The right-hand side may be further rewritten by expansion of the vector triple 
product In the first term: 

A{(ExB)*n/y 0 } * A{(v*n)(B 2 /u 0 )-(B*n_)(v_-B)/p 0 } + p(v*n)A(i?> (8) 


where the last term on the right represents the resistive contribution in 

Eq. (7). The quantity R, vanishes at station 1 where j*0, i.e., R =0. It 

1 1 
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also vanishes at station 2, i? 2 *0, except when this point Is located In the 
magnetopause Itself where j. is large. It Is difficult to produce a reliable 
theoretical estimate of i^but, for illustrative purposes, we may assume a 
scalar resistivity In which case the resistive term reduces to 


p(v/ji)Mfl} * A{-qn-VB 2 /2y 2 } 

(9) 

p(v«n)i? 2 = -r\ 2 n.*vB 2 /2y 2 

(10) 


Since the magnetic pressure usually has a minimum in the magnetopause (and 
since (v/n.)<0), we expect # 2 > 0 In the outer half of the magnetopause struc- 
ture and z? 2 <0 in the Inner half. 

If Eq. (1) is divided by the (negative and constant) mass flux p ( 
and the expressions (2), (3), (5), and (8) are used for the individual terms, 
we obtain 


4 {r} + 4 {! + "z + Q * 4( 1+ t)v } + 4 1 ^ 7^7 110 

© © © © 

In term ( v *n_) may be replaced by the Alfven speed ±(B^n)»''(1-a 1 l/u^ as 
is appropriate for a rotational discontinuity. The upper sign applies north 
of the reconnection line in the standard reconnection model where B/n<0, and 
the lower sign south of that line where B*n>0. 

In summary, for isotropic pressure, a=0, terms (]) and Q) in Eq. (11) 
represent the Increase in kinetic energy and enthalpy per unit mass of the 
plasma as It moves from station 1 to station 2 while terms ® and ® rep- 
resent the nonresistive part of the electromagnetic energy input. For 
nonisotropic pressure, a^O, some of the terms contributing to the enthalpy 
appear on the right-hand side of Eq. (11). Also, the term Q is the net heat 
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flow and is expected to be positive or zero; W 2 is the viscous contribution 
at station 2 and Is expected to be negative 1r the magnetopause and approxi- 
mately zero In the boundary layer plasma; R z Is the resistive contribution 
and Is expected to be positive in the outer half of the magnetopause, negative 
in the inner half, and zero in the boundary layer. 

Except as noted below, the terms in Eq. (11) are directly 

measured on ISEE 1 and 2 and will form the basis for our check of the energy 
balance. Apart from small viscous corrections, the perpendicular pressure, 
p A , is usually reasonably well determined by the MPE/LANl fast plasma analyzer 
throughout the magnetopausa a -id boundary layer. On the other hand, the parallel 
pressure, p^ , Is well df-rrmined only on the magnetosheath side of the mag- 
netopause where the magnetic field usually lies near the spacecraft equatorial 
plane. On the manneto spheric side, the magnetic field orientation falls in 
the blind direction of the instrument along the spacecraft spin axis. Thus 
o^, but not a 2 . Is measured directly. However, a 2 can be Inferred from the 
measured densities by use of the relation (1 -Oj ) = P 2 (^“ a 2 ^ applies 

across a rotational discontinuity. This procedure was also employed In 
Paper II. It is not strictly justified when station 2 lies in the magneto- 
pause. For this and other reasons it may occasionally yield unphysical 
results: If the inferred a 2 value is large and negative while at the same 

time the measured ratio Pj^Uq/B* is small then the formula a £ = (p /| 2 ‘ p l2^ u 0 ^ B 
will lead to a negative value for pu £ . Such points have been eliminated. 


IN <s| 
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3. OBSERVATIONS 

September 8, 1978 : Detailed information concerning this high-speed 

flow event may be found in Papers I and II. For purposes of review, the 
first three panels of Fig. 1 show the plasma density, the GSM Z component 
of the magnetic field, and the field magnitude, measured by ISEE 1 during 
its outward passage through the magnetopause region. The principal magneto- 
pause crossing is indicated by the reversal of B z at 00.44 UT. It is followed 
by a partial reentry of the spacecraft in the interval 00.46 - 00.51 UT which 
resulted from an outward -inward motion of the magnetopause. The bottom panels 
of Fig. 1 show the kinetic energy v 2 /2 and the enthalpy (5/2)p/p ■ ( 5/2) kT/m 
where p =(p |( +2p i )/3. The high-speed flow regions in the magnetopause and 
boundary layer are clearly seen. It is also observed that the high-speed 
plasma was much hotter than the ambient plasma. On the average, the enthalpy 
increase was comparable to, but somewhat smaller than, the kinetic energy 
increase. The ISEE 2 data for this event are not shown but they have the 
same qualitative properties as the ISEE 1 data and they will be used below. 

As discussed in Paper II, ISEE 2 experienced only a single crossing of the 
magnetopause during part of which a data gap occurred. 

In order to check the validity of the energy equation (11) we now 
apply that equation between the reference point 1 in the magnetosheath 
adjacent to the magnetopause and those individual measurements at a point 
2 in or inside the magnetopause which meet the selection criteria described 
in detail in the appendix of Paper II. The reference values are averages 
over several measurements. 

Figure 2 shows a scatter diagram of the terms a{{T)+( 2)} versus 
a{ / T)+ : ' 4)}, i.e., approximately, of the kinetic energy and enthalpy increase 
versus the electromagnetic energy supplied. In order to illustrate that the 
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effect of the a terms is relatively small, we have shown terms 'T) and Q 
both with the measured aj values and with <*j =0. In both cases a 2 is deter- 
mined from the density ratio p x /p 2 as described earlier. 

If the ideal reconnection model applies, then, in the absence of Q, w 2 , 
i? 2 , and experimental uncertainties, the data in Fig. 2 should fall along the 
45° line. The figure shows that, the data from the main ISEE 1 crossing and 
from the ISEE 2 crossing fall mainly below the line whereas the data from 

the partial ISEE 1 crossing fall above it. This effect may be systematic 

and, if so, we believe it can be explainer qualitatively by our neglect of 

the terms Q and as described below. 

Data from the main ISEE 1 crossing and from the ISEE 2 crossing have 
station 2 located mainly in the boundary layer or in the innermost part of 
the magnetopause. As pointed out in Section 2, for this situation one 
expects that the net heat flow Q away from the magnetopause could be sub- 
stantial (Q>0) while the viscous work w 2 and the resistive term r z are both 
small. The expected net result is that the data points should fall below 
the 45° line ir Fig. 2 as indeed they do. A quantitative estimate of Q, based 
on observed fluxes of reflected magnetosheath ions, will be given in the dis- 
cussion section. 

Data from the partial ISEE 1 crossing have station 2 located mostly in 
the outer half of the magnetopause. In this case, the net heat flow Q may 
be small and the negative term may provide the largest correction. Neg- 
lect of this term causes the data points in Fig. 2 to fall above the 45° 
line as indeed they do. In the discussion section, we shall estimate w 2 
and r 2 and show that W 2 may be substantial while R 2 (which would move the 
data points in Fig. 2 to the left) is probably small. 


In Fig. 2, the sign chosen for term (J) corresponds to a crossing north 
of the reconnection line since that was the location deduced in Papers I 
and II. Reversal of the sign of this term moves all points in Fig. 2 to 
the left as shown in Fig. 3. Such a shift would destroy the rough overall 
agreement of the data with the 45° line. Thus the energy balance supports 
the identification of this case as a crossing north of the reconnection line. 

August 9, 1978 : The format of the data presentation is the same as in 
the previous case. As shown in Fig. 4, the data used for this event cor- 
respond to a partial crossing of the magnetopause during the interval 
19.34 - 19.53 UT. During the rapid main crossing, at 20.10 UT no plasma 
acceleration was measured. It appears that an abrupt change in the magneto- 
sheath conditions may have caused the density change at 19.53 UT at the end 
of the partial crossing. For this reason, the data for the reference point 1 
was formed as an average over magnetosheath measurements irmediately pre- 
ceding the event. The increases in plasma kinetic 

energy and enthalpy were smaller than in the September 8, 1978, event and, 
on the average, the latter increase somewhat exceeded the former. 

Figure 5 shows the terms a{ ©+©> versus M(3)+(4)} for the August 9 
event (for ISEE 1 and 2). All of the data points are seen to fall above the 
45° line. This behavior is consistent with our expectation for data from 
partial magnetopause crossings where static 2 is located in the magneto- 
pause so that ^ 2 <0. If this term were included, the data points in Fig. 4 
would move down towards the 45° line. For this event there was no clear 
evidence of reflected magnetosheath ions so that the ion heat flow away 
from the magnetopause was presumably small. And the R 2 term should be 

small because vB 2 is small. 

— 2 
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In Fig. 5 the sign of term © chosen corresponds to a crossing south 
of the reconnection line since that was the location deduced in Paper II. 
The opposite choice of sign is made In Fig. 6, in which most of the points 
are seen to have negative values of the abscissa. It would be difficult 
to argue that the data in this figure have an tendency to agree with the 
45° line. Thus the energy balance supports the identification of this 
case as a crossing south of the reconnection line. 
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4. DISCUSSION 

Figure 2 shows that the data points from the complete magnetopause 
crossings by ISEE 1 and ISEE 2 on September 8, 1978, fall systematically 
below the 45° line. In other words, these data Indicate that ail of the 
electromagnetic energy dissipated In the magnetopause does not show up 
as kinetic energy and enthalpy Imparted to the plasma as It crosses this 
current layer. We shall now show that most of the missing energy leaves 
the system In the form of heat flow associated with a population of magneto- 
sheath ions which Is reflected and energized at the magnetopause and returns 
upstream as a beam flowing nearly along the magnetic field. As reported In 
Paper II, this beam has a density of approximately 20% of the magnetosheath 
density and a velocity of about 500 km/s. It Is easy to show that the heat 
flow vector In a plasma consisting of two cold Ion beams, of density rp 
and (1-rJpj, having velocity difference Av, Is 

x 1 «> 1 r(l-r)(l-2r)|Av| 2 Av 

where In the present case Av * - |av |B / |B 1 | . From Eq. (5) with 3^=0 and 
(y_ x -n.) * (B.n)/rr^J7v7 we then find 

Q ' ‘ * ‘s r t 1 - r )0-2 r )|av| 3 /v A i 

where v^ * l-i I ^ " a i )/u 0 p i * For r * °-2> I^Y-I * 500 km/s and v Al * 300 km/s 
we find Q ■ 2xl0 14 cm 2 /s 2 . Values of Q of this order of magnitude are also 
obtained by direct integration over the measured Ion distribution function. 

If the data points for the two complete magnetopause crossings In Fig. 2 
are moveH upwards by 2*10 14 cm 2 /s 2 then, on the average, their location will 
coincide well with the 45° line. The remaining discrepancy may be easily 
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accounted for In terms of the uncertainties in the measurements. 

It Is more difficult to provide reliable estimates of the terms W 2 
and # 2 which, as we have argued, tend to move data points In Figs. 2 and 
4 from the partial crossings downwards and to the left. However, we may 
write these terms as 

w * D 

2 v h(v_*n_) 

(AV .) 2 

R z * D n h(vnj 

where h Is the magnetopause thickness and av 2 and avJ are characteristic 
changes of flow speed and Alfvdn speed In the magnetopause. Also 0 y and D n 
are the viscous (m/p) and resistive (n/u Q ) diffusion coefficients, respectively. 
As an example, for the September 8 event, we may choose Av * 400 km/s, 

Av^ « 150 km/s, h * 100 km, and (y_*r[) * 30 km/s. In order to produce a 
value of 3x1 0 14 cm 2 /s 2 for w 2 or J? 2 , which is of the order of the corrections 
needed (see Fig. 2), the diffusion coefficients would have to be D v -5xl0 12 cm 2 /s 
and D^-4*l 0 13 cm 2 /s . These values should be compared to the diffusion coeffi- 
cient D ■ R 2 /T corresponding to a diffusion speed of one gyro radius R per 
Ion gyro period T . For R “50 km and T -is we have D * 2.5xlO l3 cm 2 /s. Comparison 
with D now Indicates that the above value for D y is reasonable while that for 
D r is too large. Thus we expect the dominant correction term to be With 
reasonable assumptions, this term can become sufficiently large to move the 
data points for the partial crossings down to the 45° line In Fig. 2. 

For the August 9 event, the quantity (y^-n.) Is not known. Therefore It 
becomes Impossible to produce reliable estimates of the required diffusion 
coefficients. However the magnetic field magnitude change was small during 
the crossing so that the correction term R z was probably small In this case also. 
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In summary, our examination of the energy balance during two plasma 
acceleration events (September 8 and August 9, 1 978) has led to the fol- 
lowing results: 

(1) To within experimental and other uncertainties the data are 
consistent with the Interpretation in terms of magnetic field 
reconnection offered In Papers I and II. This consistency 
Includes an explicit check on the sign of the magnetic 
field component normal to the magnetopause. 

(11) The enthalpy Increase of the magnetosheath plasma as It 

crosses the magnetopause during reconnection Is substantial 
and is comparable to Its Increase In kinetic energy. 

(Ill) During the September 8 event, a substantial Ion heat flow 
from the magnetopause outward Into the magnetosheath was 
present. This heat flow carries away some 20% of the 
total electromagnetic energy dissipated In the magnetopause. 

(1v) In the magnetopause structure Itself, the energy budget may 
Include substantial contributions from viscous effects. 

The experimental uncertainties Involved In this examination of the 
energy budget are large. And, as emphasized In Paper II, the theoretical 
model to which the data are compared Is undoubtedly oversimplified In many 
respects. Nevertheless, In our view, the results reported here add sub- 
stantial support to the Interpretation of these two plasma acceleration 
events In terms of the magnetic field reconnection process. We have also 
established that the energy conversion In the magnetopause during reconnection 
Is not purely electromechanical. A substantial amount of dissipation occurs 
and produces heating of the plasma as it crosses the current layer. Such 
behavior, while perhaps unexpected, is nevertheless entirely consistent with 


15 


the jump conditions across rotational discontinuities In a plasma with 
nonisotropic pressure (Hudson, 1970; 1971; 1973). 
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FIGURE CAPTIONS 

Fig. 1 Plasma and magnetic field data from ISEE 1 for a 36-min. interval 
on September 8, 1978, From the top, the panels show: the plasma 

density N p ; the GSM Z component of the magnetic field the field 
magnitude B; the kinetic energy v 2 /2 and the enthalpy 5kT/2m, per 
unit mass (both in units of 10 14 cm 2 /s 2 ) . Note the high values of 
the latter two quantities in the magnetopause (MP) and boundary layer 
(BL) regions. Heating in the magnetosheath adjacent to the magneto- 
pause is also evident. 

Fig. 2 Scatter plot of terms (T) and (2) versus terms © and (J) in the 
energy balance equation (11) for the September 8, 1978, event. The 
sign of term © corresponds to a crossing north of the reconnection 
line (B n <0). The symbols (+) and (•) denote the principal magneto- 
pause crossing by ISEE 1 and the only crossing by ISEE 2, respectively. 
The symbols (X) and (0) denote the left-hand and right-hand halves 
of the partial ISEE 1 crossing. Bars extending from these symbols 
indicate the change in the location of data points if one puts <^=0. 

Fig. 3 Same as Fig. 2 but with the opposite sign of term (4). 

Fig. 4 Plasma and magnetic field data from ISEE 1 for a 60-min. interval 
on August 9, 1978. See Fig. 1. 

Fig. 5 Scatter plot of terms (T) and © versus terms (T) and (£' in the 
energy balance equation (11) for the August 9, 1978, event. The 
sign of term (4 corresponds to a crossing south of the reconnection 
line (B n >0). 

Fig. 6 Same as Fig. 5 but with the opposite sign of term (A 
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An Analytic Model of Current Layers in a Collision-Free 
Plasma With Anisotropic Pressure 

Keith B. Kirkland 
M.S. Thesis 
September 1980 

Abstract 
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In any ionized gas, or plasma, with an imbedded magnetic field, an abrupt 
change in the orientation and strength of the field at a surface implies the 
existence of a thin sheet of current at that location. Such current sheets 
provide boundaries between different magnetic domain and are cornnon in the 
solar system. For example, they are found in solar flares and in planetary 
"bow shocks" and magnetopauses (which occur where the solar wind blows magnetic 
field lines and plasma of solar origin past planetary magnetic fields). Further 
more, under certain conditions, ions moving in the current sheets may be acceler 
ated to very high energies. Such acceleration occurs sporadically during 
solar flares and geomagnetic storms. 

Current sheets may be divided into four categories: shocks, tangential, 

rotational, and contact discontinuities. Analytic solutions describ.ng the 
detailed self-consistent behavior of the magnetic field and the charged par- 
ticles have been sought for current sheets with a nonvanishing normal magnetic- 
field component in a collisionless plasma. Tangential dis^v itinuities, which 
have no normal components, have not been dealt with. In the cold plasma limit, 
analytic models to date have assumed the plasma to consist of a single beam 
of cold ions and an accompanying cold electron beam, implying vanishing plasma 
pressure. This leads to solutions for shock solitons (for which, conditions 
are identical on the two sides of the layer) but not to rotational or contact 
discontinuities. 

In this thesis, in order to include finite pressure effects, two cold ion 
beams with different energy were used. The resulting plasma has a finite pres- 
sure along the magnetic field but no pressure across it. With this modifica- 
tion, solutions for contact discontinuities were found and rotational dis- 
continuities appear to be possible as well as shock solitons. Particular 
attention was paid to a class of contact discontinuities for which the mag- 
netic field is coplanar, i.e., it is confined to a plane normal to the current 
layer. Among these solutions, four subclasses were found: (i) soliton-like 

layers in which the magnetic field and plasma are identical on the two sides; 

(ii) periodic solutions consisting of an infinite array of current layers; 

(iii) reflection layers, in which the incident ion beam is reflected; ( i v ) 
trapped layers, where the ions are trapped within the sheet. Finally, the 
ground has been prepared for exploration of contact discontinuities where 
B is not coplanar, rotational discontinuities and shock solitons. 


A Stability Analysis of the Interface Between the Magnetopause Boundary 
Layer and the Magnetosphere Including Coupling to the Ionosphere 

Susan E. Minas 
M.S. Thesis 
September 1981 

Abstract 


The sun is constantly emitting a flux of charged particles, or plasma, 
known as the solar wind, which flows out into the solar system and past the 
earth. The magnetic field of the earth carves out a cavity, referred to as 
the magnetosphere in this flow. Some of the flowing plasma makes it across 
the cavity boundary, called the magnetopause, and continues to flow in the 
general direction of the solar wind, in a thin boundary layer. 

The plasma flowing in the boundary layer moves across the magnetic field 
lines of the earth, setting up a natural MHD generator. It is convenient to 
consider this generator to be hooked up to the resistive auroral ionosphere, 
with the magnetic field lines acting like wires. Phenomena occurring in the 
boundary layer are expected to have noticeable effects in the auroral region. 

It is known from fluid mechanics that the interface between a flowing 
and a stationary fluid is susceptible to the Kelvin-Helmhol tz (K-H) instabi- 
lity. This is true for a plasma interface as well with the additional feature 
that any shear in the magnetic field at the interface wends to suppress the 
instabil ity. 

A region of the low latitude boundary layer (LLBL), on the morning side 
of the earth, within a distance H from the equatorial plane is modeled as an 
MHD generator for the ionospheric load. In the case of perfect coupling to 
the ionosphere, magnetic field lines are assumed to be equi potentials and 
any potential developed across the boundary layer is impressed on the iono- 
sphere, causing a current to flow there which becomes part of a current loop. 
This current flows across the LLBL and into and out of both the LLBL and 
ionosphere along magnetic field lines. Perturbations in the field aligned 
current due to electric field perturbations felt in the ionosphere are 
matched to the perturbation current in the LLBL associated with perturba- 
tions in the magnetic field there. All perturbations are expressed as series 
expansions in Z, the distance from the equatorial plane, and the dispersion 
relation is arrived at keeping terns to lowest significant order in Z. 

The dispersion relation to this order indicates that the coupling to 
the ionosphere decreases the growth rate of the instability. It does not 
exhibit any threshold in flow speed below which the interface is stable. 

Using obstr/ed values for the model parameters it is shown that the growth 
rate remains sufficiently large so that large amplitude waves would develop 
during the time required for the K-H waves to travel from the subsolar point 
to an observation site in the mor.iingside boundary layer. Thus, it seems 
reasonable to believe that the Kelvin-Helmhol tz instability at this inter- 
face is responsible for the observed modulation in boundary layer thickness. 

A brief summary is given of the steps needed to carry the Z expansion 
to the next order. 


Spectral Analysis of Magnetic Fields Around 
the Earth's Magnetopause 

Fernando Fonseca 
M.S. Thesis 
December 1983 


Abstract 


Spectral analysis was performed on magnetic field data obtained near 
the earth's magnetopause by the fluxgate magnetometers onboard the satellites 
in the International Sun-Earth Explorer mission (ISEE). Data were analyzed 
from the plasma boundary layer just inside the magnetopause as well as from 
the adjoining regions of the magnetosphere and the magnetosheath. A survey 
of the techniques available to analyze many short segments of data was 
performed and the most appropriate methods were selected and Implemented. 

The analysis was performed in the frequency range .1 to 1 hertz. Variations 
in the power spectra of the magnetic field fluctuations of the form f‘ 2 
were found in the magnetosheath and in the boundary layer. Study of the 
spectral matrices indicated that the noise in these regions has a propagation 
direction approximately along the magnetic field and may be associated with 
the whistler mode. In the magnetosphere a power spectrum of the form f -3 
was obtained. Except for a small peak at the satellite spin frequency 
observed in the magnetosphere, no signficant peaks were found in the power 
spectra from any of the three regions. 


